Nonpositive Curvature: a Geometrical Approach to 
Hilbert-Schmidt Operators* 

Gabriel Larotonda^ 

Abstract 

We give a Riemannian structure to the set E of positive invertible unitized Hilbert-Schmidt operators, 
by means of the trace inner product. This metric makes of E a nonpositively curved, simply connected 
and metrically complete Hilbert manifold. The manifold S is a universal model for symmetric spaces 
of the noncompact type: any such space can be isometrically embedded into E. We give an intrinsic 
algebraic characterization of convex closed submanifolds M. We study the group of isometries of 
such submanifolds: we prove that Gmi the Banach-Lie group generated by M, acts isometrically and 
transitively on M. Moreover, Gm admits a polar decomposition relative to M, namely Gjv/ ~ M X if as 
Hilbert manifolds (here K is the isotropy of p = 1 for the action /g : p i— > gpg*), and also Gm/K ~ M 
so M is an homogeneous space. We obtain several decomposition theorems by means of geodesically 
convex submanifolds M . These decompositions are obtained via a nonlinear but analytic orthogonal 
projection Hm '■ — > Ai, a map which is a contraction for the geodesic distance. As a byproduct, we 
prove the isomorphism NM ~ E (here NM stands for the normal bundle of a convex closed submanifold 
M). Writing down the factorizations for fixed e"^, we obtain e" = e^e"e^ with e^ S M and v orthogonal 
to M at p = 1. As a corollary we obtain decompositions for the full group of invertible elements 
G ~ M X cxp(TiM-L) X i^Q 



1 Introduction 

The aim of this paper is to relate the algebraic and spectral properties of the Banach algebra of 
unitized Hilbert-Schmidt operators, with the metric and geometrical properties of an underlying 
manifold E. This is a paper on applied nonpositively curved geometry because we first show how 
the familiar properties of the operator algebra translate into geometrical notions, and then we use 
the tools of geometry in order to prove new results concerning the operator algebra. 

In this paper we study the cone of positive invertible Hilbert-Schmidt operators (extended 
by the scalar operators) on a separable infinite dimensional Hilbert space H . The metric in the 
tangent space at the identity is given by the trace of the algebra. The local structure induced 
by the metric is smooth and quadratic; it can be situated in the context of the theory of infinite 
dimensional Riemann-Hilbert manifolds of nonpositive curvature (c/. Cartan-Hadamard manifolds, 
as introduced by Lang [T^, McAlpin [22 , Grossman jl^ and others). It is then a paper on 
Riemannian geometry. On the other hand, since the manifold S is clearly not locally compact, 
some of the standard results for Hadamard manifolds require a different approach. The geometry 
is then related to the geometry of the metric spaces in the sense of Aleksandrov [5] . It turns out 

*2000 MSG. Primary 22E65; Secondary 58E50, 53C35, 53C45, 58B20. 

tinstituto de Ciencias, Universidad Nacional de General Sarmiento, JM Gutierrez 1150 (1613) Los Polvorines, 
Buenos Aires, Argentina, e-mail: glaroton@ungs.edu. ar 

^Keywords and phrases: Hilbert-Schmidt class, nonpositive curvature, Banach-Lie group, homogeneous man- 
ifold, operator decomposition 



that the notion of convexity (together with the fact that E is a simply connected and globally 
nonpositively curved geodesic length space) plays a key role in our constructions. It is then a paper 
on metric geometry. 

Through the years, several authors have studied the relationship of geometry and algebra in sets 
of positive operators, with different approaches that led to a variety of results. In his 1955's paper 
G.D. Mostow gave a Riemannian structure to the set M+ of positive invertible matrices; the 
induced metric makes of a nonpositively curved symmetric space. Mostow showed that the 
algebraic concept behind the notion of convexity is that of a Lie triple system, which is basically 
the real part of a given involutive Lie algebra g. The geometry of bounded positive operators in 
an infinite dimensional Hilbert space was studied by G. Corach, H. Porta and L. Recht |Tn]|13|P5] 
among others, using functional analysis techniques. This area of research is currently very active 
(see [8] |9] for a list of references) . 

1.1 Main results 

In this paper we study the geometry of a Hilbert manifold E which is modeled on the operator 
algebra Tic of unitized Hilbert-Schmidt operators. In Section [2] we introduce the objects involved 
and prove some elementary results. The manifold E is the set of positive invertible operators of 
He- Let He* be the classical Banach-Lie group of invertible (unitized) Hilbert-Schmidt operators 
|17| . The manifold E has a natural Tic'-invariant metric < x,y >p = < xp~^ ,p~^y which makes 
it nonpositively curved (we define < a + a,/3 + 5 >2= oi[3 + Atr{h*a) whenever a,/3 G C and a, 6 
are Hilbert-Schmidt operators). Let e^ and exp(x) stand for the usual analytic exponential, i.e. 
exp(a;) — X]n>o fj- This map is injective when restricted to TYk, the set of self-adjoint operators. 
Let ln(p) stand for its real analytic inverse. We have exp(?-^R) = E C Hk, and the exponential 
map induces a diffeomorphism onto its image, so we identify the tangent space at any point of the 
manifold S with the set of self-adjoint operators T^r, namely TpE ~ Tia for any p e E. In Section 
[3] we prove 

Theorem A: Forp,q G E, the geodesic obtained from Euler's equation by solving Dirichlet's 
problem is the smooth curve Jpq{t) — p'^ [p^^ qp^^Yp^ , hence 

/ \ 1 /—I — l\i 

ExPp(w) = p^ exp(p 2 y p 2jp2 

is the Riemannian exponential ofY., for any v £ TpE. Both Expp : TpE — > E and its differential 
map d (Expp) ^ : TpE TExpp(i))E are -isomorphisms for any p G E, w G T^E. The curve ^pq is 
the shortest piecewise smooth path joining p to q, hence 

dist{p,q) = \\ln{p^q-^pi)l, 

is the distance in E induced by the Riemannian metric. The metric space (E, dist) is complete, and 
it is globally nonpositively curved. 

The curve obtained via Calderon's method of complex interpolation fTj between the quadratic 
norms || • ||p and || • ||^ is exactly the short geodesic in E joining p to q (the proof of [1] can be 
adapted almost verbatim). 

In [15] . N. Grossman proves that the inequality 

|M(Expp)jz«)l|^>||«;L^^^,„, (1) 
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leads to the minimality of geodesies in a simply connected, complete Hilbert manifold. This ap- 
proach is also carried out by McAlpin [22j . The following operator inequality involving the differ- 
ential of the usual exponential map 

||e--/2dexpjy)e--/2||, >|l2/||, (2) 

is the translation to our context of the inequality ^ above. The convexity of Jacobi fields can be 
deduced from the non positiveness of the sectional curvature, hence the proof of eqn. ([2]) stems in our 
context from the Cauchy-Schwarz inequality for the trace inner product. We follow the exposition 
of Lang [50] on this subject. On the other hand, ^ can be proved with a direct computation 
[6]. With this approach the metric completeness of the tangent spaces is not relevant: in Theorem 
3.1 of [5], the authors prove the minimizing property of the geodesies in a non complete manifold. 
The inequality above, in our context, can be also interpreted as the Hyperbolic Cosine Law (see 
Corollary [332]) 

> 6^ + - 2feccos(Q;). 

Here a, b, c are the lenghts of the sides of any geodesic triangle in S, and a is the angle opposite to 
a. From this inequality also follows that the sum of the inner angles of any geodesic triangle in S 
is bounded by tt. 

If ^ is a set of operators, we use to denote the set of positive operators of A; note that 
("^C*)^ = 1]. In Section |4] we show that a submanifold M C S is geodesically convex if and 
only if its tangent space at the identity m is a Lie triple system. Clearly any such submanifold is 
nonpositively curved, and Theorem 14. 181 states: 

Theorem B: For any geodesically convex, closed submanifold M ~ exp(m) C S there ex- 
ists a connected Banach-Lie group Gm = (exp(m © [m,m])) C He* which acts isometrically and 
transitively on M. Moreover, the polar decomposition of the elements of Gm reduces to M in the 
sense that G^j = M . Let K be the isotropy of 1 for the action; then K is a connected Banach-Lie 
subgroup of Gm cind there is an isomorphism Gm — M x K. Ln particular any convex submanifold 
M of S is an homogeneous space for a suitable Banach-Lie group, which is an analytical subgroup 
of 'He* ■ The submanifold M is flat if and only if M = Gm is an abelian Banach-Lie subgroup of 
He' ■ 

The existence of smooth polar decompositions for the involutive Banach-Lie groups can be 
obtained from the general results of Neeb ([24], Theorem 5.1). Neeb introduces the notion of 
seminegative curvature (SNC) on Banach-Finsler manifolds M, given by the condition of inequality 
([T]) above, plus the condition that d (Exp^)^^ should be invertible for any v e TpAI (the metric of M 
sould be invariant under parallel transport along geodesies). Neeb proves (Theorem 1.10 of [24]) that 
in a connected, geodesically complete manifold with SNC, the exponential map Exp^ : TpM —>■ M 
is a covering map and M is metrically complete, a result which extends that of Grossman and 
McAlpin mentioned above to the Banach-Finsler context. 

The manifold E can be decomposed by means of any convex closed submanifold M. Let NM 
be the normal bundle of M . In Section 5 we prove 

Theorem C: For any convex closed submanifold Af C S there is a nonlinear, real analytic 
projection Hm : S — > M , which is Hm is contractive for the geodesic distance 

dist(HM(p), T^Miq)) < dist(p, q) for any p,q (ET.. 

The point UmIp) is the (unique) point of M closest to p. It can also be viewed as the unique point 
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in M such that there exists a geodesic through p orthogonal to M at Hm (p) ■ The exponential map 
{p,v) ^—^ Expp(?;) induces an analytic Riemannian isomorphism NM ~ S. 

Since UmIp) is the point in M closest to p, one can prove the existence of such a point using a 
metric argument vahd in any nonpositively curved geodesic length space |18j . We choose to give a 
differential-geometry argument here. 

In Section 6 we exhibit a decomposition for the submanifold M = A of positive diagonal 
operators, which is a maximal abelian subalgebra of Tic. This decomposition theorem (Theorem 
16. 2p takes the form of a factorization e° = de^d, where v has null diagonal and d is an invertible 
diagonal operator. We stress that there is no known algorithm that allows to compute d explicitly 
(not even if we reduce the problem to 3 x 3 matrices, that can be thought of as a particular case of 
the general theory). As a corollary to the decomposition theorems we obtain 

Theorem D: Any invertible operator g £ Tic* admits a unique polar decomposition relative 
to a fixed closed convex submanifold M = exp(m). Namely g — e^e"u where a; £ m, w G and 
u e U{Ti.c) is a unitary operator. The map g i—>- {e^,e'",u) is an analytic bisection which gives the 
isomorphism 

He' ^ M X exp(m-L) x UiHc)- 

This isomorphism generalizes the decomposition of given in [llj . 

In Section 7 we show that the manifold E can be decomposed by means of a foliation {Ea}a>o 
of totally geodesic submanifolds, namely 

I]= UE\= U{a + AGE, a = a*a Hilbcrt Schmidt operator). 

A>0 A>0 

There is a Riemannian isomorphism E ~ Ei x R>o induced by the projection Tlsi of Theorem C 
above. As an application, we show a decompositon relative to the algebra of positive invertible 
nx n matrices: fix an n-dimensional subspace S C H, let be the orthogonal projection to S and 
— 1 ~ Pg the orthogonal projection to S*^. Let B{S) stand for the algebra of bounded linear 
operators of S. Let R e B(S')+ ~ M^, and consider the set 



= 



B{S^) a Hilbert-Schmidt operator , Y e B(5, 



Let U{Ti.c) be the Banach-Lie subgroup of unitary operators in Tic* ■ 

Theorem E: For any g G Ti-c' there is a unique factorization g = Xre'"u where A G K>0; 
u G U{Ti.c) is a unitary operator, r ~ RPs + Qs o*^^ v £ V. In particular 

He' ^ A/+ x exp(o) X R>o x W(Hc). 

The manifold E can be regarded as a universal model for the symmetric spaces of the noncompact 
type, namely 

Theorem F: For any finite dimensional real symmetric manifold M of the noncompact type 
(i.e. with no Euclidean de Rham factor, simply connected and with nonpositive sectional curvature) , 
there is an embedding M ^ Ti which is a diffeomorphism between M and a closed geodesically convex 
submanifold of T,. If we pull back the inner product on E to M , this inner product is a positive 
constant multiple of the inner product of M on each irreducible de Rham factor. 



4 



The proof of the theorem is straightforward fixing an orthonormal basis of H (see Section 17. ip 
and recaUing the well known result [14 that for any such space M there is an almost isometric 
embedding of M into GL{q)^ , where g is the Lie algebra of the Lie group Iq{M) (the connected 
component of the identity of the group of isometrics of M) . 

2 Background and definitions 

Let B(_ff) be the set of bounded operators acting on a complex, infinite dimensional and separable 
Hilbert space iJ, and let HS be the bilateral ideal of Hilbert-Schmidt operators of B(i?). Recall that 
HS is a Banach algebra (without unit) when given the norm ||a||2 = tr{a*a)^ (see [57] for a detailed 
exposition on trace-class ideals). We will use HS'' to denote the closed subspace of self-adjoint 
Hilbert-Schmidt operators. In B(i7) we define 

nc = {a + X: a e HS, A G C}, 

the complex linear subalgebra consisting of Hilbert-Schmidt perturbations of scalar multiples of the 
identity (the closure of this algebra in the operator norm is the set of compact perturbations of 
scalar multiples of the identity) . There is a natural Hilbert space structure for this subspace (where 
scalar operators are orthogonal to Hilbert-Schmidt operators) which is given by the inner product 

<a + \b + P Atr{ab*) + Xp. 

The algebra Tic is complete with this norm. The model space that we are interested in is the real 
part of He, 

nR = {a + X: a* a, a e HS, A G M}, 

which inherits the structure of (real) Banach space, and with the same inner product, becomes a 
real Hilbert space. 

Remark 2.1. By virtue of trace properties, < xy,y*x* >^=< yx,x*y* >^ for any x,y G Tic, and 
also < zx, yz >2—< xz, zy >2 for x,y & He and z G Hs.. 

Let E := {A > : A E Hm} be the subset of positive invertible operators in Hr. It is clear that 
E is an open set of T^r (for instance, using the lower semi continuity of the spectrum). 

Remark 2.2. For p G E, we identify TpE with Hm, and endow this manifold with a (real) Rieman- 
nian metric by means of the formula 

<x,y >p-=< p~^x,yp~^ >2=< xp~^,p~^y >2 . 
Throughout, let ||a;||p :=< x,x >J. Equivalently, \\x\\^ = \\p^^ xp^^W^. 

Lemma 2.3. The covariant derivative in E (for the metric introduced in Remark \2.2\) is given by 
{Vx^lp - {X{Y)}p - i {Xp p-' Yp + Yp p-' Xp) . (3) 

Here X{Y) denotes derivation of the vector field Y in the direction of X performed in the linear 
space Hr. 

Proof. Note that V is clearly symmetric and verifies all the formal identities of a connection; 
the proof that it is the Levi-Civita connection relays on the compatibility condition between the 
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connection and the metric, < X,Y >^=< VjX,Y >^ + < X, V^F >^ (see for instance pp] 
Chapter VIII, Theorem 4.1). Here 7 is a smooth curve in E and X,Y are tangent vector fields 
along 7. This identity is straightforward from the definitions and the properties of the trace. □ 

Let r" = e"'"('^) (here r G E, a G M). The exponential is given by the usual series; note that any 
positive invertible operator has a real analytic logarithm, which is the inverse of the exponential 
in the Banach algebra. Note that aba > whenever a,b > and also r" > whenever r > and 
a G R. 

Euler's equation V^7 = for the covariant derivative introduced above reads 7 = 77 7, and 
it is not hard to see that the (unique) solution of this equation with 7(0) = p, 7(1) ^ q is given by 
the smooth curve 

Ipq {t)=P^ (P" hp'^fp^- (4) 

Remark 2.4. We will use Exp^ : TpS ^ E to denote the exponential map of E. Differentiating at 
t = the curve above, we obtain jpq{0) ~ \n{p^ i qp^ i )pi , hence 

—1/ \ -1 / — - — -\ - / \ - / — - — -\ - 

Expp (q) = p2 \n[p ^ qp ^)p^ and F,xpp[v) — p^ exp(p ^vp ^)p^. 

Note that by the construction above the map Expp : TpE — > E is surjective (for given q G E 
take V — pi \n{p^i qp^i)pi , then Expp(w) = q). Rearranging the exponential series we get the 
expressions Expp(u) — p — e^^ ^p. 

Lemma 2.5. The metric in E is invariant under the action of the group of invertible elements: if 
g is an invertible operator in Tic, then Ig{p) — gpg* is an isometry ofT,. 

Proof. First note that for any ip ^ H we have < gpg*'ip,'4> >=< pg*"!!), g*'il' >=< pr],ri > > 
assuming p > and g invertible, so Ig maps E into itself. Also note that d{Ig)r(x) = gxg* for any 
X G T^E, hence 

Il5a;.9*ll',,. -^<gxg*{g*)-^r-^g-^,{g*)-\-^g-^gxg* >,= 

= < gxr-^g-\ {g*)-^r'^xg* >,=< xr-\r~^x >,= \\x\\l 
where the third equality in the above equation follows from Remark 1 2. II □ 



3 Local and global structure 
3.1 Curvature 

We start showing that curvature in this manifold is a measure of noncommutativity, and then give 
a few definitions, which are necessary because of the infinite dimensional setting. Let [ , ] stand for 
the usual commutator of operators, [x, y] ~ xy — yx. 

Proposition 3.1. The curvature tensor for the manifold E is given by: 

y\p{x,y)z = - jP [[p^^x,p^'^y] ,p^^z] . (5) 

Proof. This follows from the usual definition £H(x, y) — V^Vj, — V^Va; — '^[x.y]- The formula for V 
given in Lemma 12.31 □ 
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Definition 3.2. A Riemannian submanifold M <Z T, is flat at p ^ M if the sectional curvature 
vanishes for any 2-subspace of TpM . The manifold M is flat if it is flat at any p € M . The 
manifold M is geodesic at p G M if geodesies of the ambient space starting at p with initial velocity 
in TpM are also geodesies of M . The manifold M is a totally geodesic manifold if it is geodesic at 
any p G M. Equivalently, M is totally geodesic if any geodesic of AI is also a geodesic ofT,. 

Proposition 3.3. The manifold S has nonpositive sectional curvature. 

Proof. Let x,y ^ TpY^. Let x — p^^xp~^ , y — p^^yp~^ . We may assume that x, y are orthonormal 
at p. A straightforward computation shows that 

Sp{x,y) ^p{x,y)y,x {< xy'^,x >^ -2 < yxy,x >^ + < y'^x,x>^} . 

Since x,y & Ti-R, x = X + a and y = f3 + b for A, /3 G K. and a, 6 G HS. The equation reduces to 

Spix, y) = -\ {tr{aH^) - tr{{ab)')} . (6) 

Note that [z,w] = tr(w*z) is an inner product on HS, so we have the Cauchy-Schwarz inequality 
tr(w*z) < tr^(w*w) tr^(z*z). Putting w = ba, z = ab, we obtain 

tr{{abf) = tr{abab) = tr{{ba)*ab) < tr^ (abba) tri (baab) = tr{a'^b'^). □ 

Proposition 3.4. Let M (ZT, be a submanifold. Assume that M is fiat and geodesic at p M . If 
x,y d TpM, thenp~"2xp^^ commutes with p^ ^ yp^ ^ . 

Proof. Since M is geodesic at p, the curvature tensor is the restriction of the curvature tensor of S, 
so in equation ^ above the right hand term must be zero if AI is flat at p. But the Cauchy-Schwarz 
inequality is an equality only if the vectors arc linearly dependent; in the notation of the previous 
theorem, we have ab — z = avu — aba for some a G M; replacing this in the above equation we 
obtain a = 1, namely ab = ba. Recalling the deflnitions for a and b we obtain the assertion. □ 

3.2 Convexity of Jacobi fields 

Let J{t) be a Jacobi field along a geodesic 7 of E, i.e. J is a solution of the differential equation 

i??j + 5n^(J,7)7 = (7) 

where Dt — is the covariant derivative along 7. We may assume that J(t) is non vanishing, 
hence 

|lJ|i3 ^< J,J>|^ - < DtJ,J + < J,J >^ < DtJ,DtJ 
- <J,J>^ < fR^(J,7)7, J >^ . 

The third term is clearly positive and the first two terms add up to a nonnegative number by the 
Cauchy-Schwarz inequality: < DtJ, J >^ < < DtJ, DtJ >^ < J, J . In other words, the smooth 
function t ^< J,J>^= \\J\\^ is convex, exactly as in the finite dimensional setting. 

3.3 Tlie exponential map 

We present two theorems that, in this infinite dimensional setting, stem from McAlpin's PhD. Thesis 
(for a proof see [22] or Theorem 3.7 of Chapter IX in [2^). First, if one identifies the Riemannian 
exponential with a suitable Jacobi lift, one obtains 
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Theorem 3.5. The map Expp : TpS — > E has an expansive differential: 

This result implies that the differential of the exponential map is injective and has closed range. 
Playing with the Hilbert structure of the tangent bundle and using the well known identity for 
operators Ker{A)-^ ~ Ran{A*), it can be proved that this map is surjective, moreover 

Corollary 3.6. The differential of the Riemannian exponential (i(Expp)„ : TpS Te^p is a 
linear isomorphism for any v £ TpT,. Hence, Expp : TpM T, is a -diffeomorphism. 

The last assertion is due to the fact that the map Exp^ : TpE ^ S is a bijection (see Remark 
above). 

3.4 The shortest path and the geodesic distance 

The following inequality is the key to the proof of the fact that geodesies are minimizing. It was 
proved by R. Bhatia [S] for matrices, and his proof can be translated almost verbatim to the context 
of operator algebras with a trace, see [1] . However since the Riemannian metric in S is complete, 
the inequality can be easily deduced from the fact that the norm of a Jacobi field is a convex map 
(in Theorem 13.51 put p = 1, v = x and w = y): 

Corollary 3.7. If d exp^ denotes the differential at x of the usual exponential map, then for any 
x,y <eHr 

IMexp^(y)lle- = l|e"*o'exp^(2/)e"^ll2 > WvL- 
As usual, one measures length of curves in E using the norms in each tangent space, 

L{a)^ r \\a{t)\l^^^dt. (8) 
Jo 

We define the distance between two points p, (7 G S as the infimum of the lengths of piecewise 
smooth curves in S joining p to q, 

dist(p, q) = inf {L{a) : q C S, a{0) — p, a{l) — q} . 

Recall (Remark 12 .41 and the paragraph above it) that for any pair of elements p, g G E, we have the 
smooth curve jpq C S, jpq(t) = pi {p~i qp~iypi joining p to q, which is the unique solution of 
Euler's equation in S. Computing the derivative, we get 

The minimality of these (unique) geodesies joining two points can be deduced from general consid- 
erations [15] ■ we present here a direct proof. 

Theorem 3.8. Let p,q G Y,. Then the geodesic jpq is the shortest curve joining p and q in S, if 
the length of curves is measured with the metric defined above 

Proof. Let a be a smooth curve in E with a{Q) — p and a{l) = q. We must compare the length of 
a with the length of jpq. Since the invertible group acts isometrically, it preserves the lengths of 
curves. Thus we may act with p^i , and suppose that both curves start at 1, or equivalently that 
p = 1. Therefore Jiq{t) :— j{t) = e*^, with x ~ hiq. The length of 7 is then ||a;||2- The proof 
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follows easily from the inequality of Corollary 13.71 Indeed, since a is a smooth curve in E, it can 
be written as a (t) = e''^*), with P = In a. Then f} is a smooth curve of self-adjoint operators with 
f3{0) = and (3{1) = x. Moreover, 

L{l) = = " Oil, = II m dt\l < f 11/3(011, dt. 

Jo Jo 

On the other hand, by the mentioned inequality, 

< lle-"^ dexp^(,)(/3(0)e-^|l, = l|dexp^(,)(/3(t))||_^^,,, = ll"WIL(„- □ 

Remark 3.9. The geodesic distance induced by the metric is given by 

dist(p,g) = ||ln(p5 q-^pi)\\,_. 

Hence the unique geodesic joining p to g is also the shortest path joining p to q. This means that 
(S,dist) is a (not locally compact) geodesic length space in the sense of Aleksandrov and Gromov 
[5]. These curves look formally equal to the geodesies between positive definite n x n matrices, 
when this space is regarded as a symmetric space. 

Corollary 3.10. //7, S are geodesies, the map / : R — > M, i dist(7(i), <5(t)) is convex. 

Proof. The distance between the points ^{t) and 5{t) is given by the geodesic at{s), which is 
obtained as the s variable ranges in a geodesic square h{s, t) with vertices {j{to),S{to), 7(^1), ^(^i)} 
(the starting and ending points of 7 and 6). Taking the partial derivative along the direction of s 
gives a Jacobi field J(s, t) along the geodesic Ps{t) — h{s, t) and it also gives the speed of at- Hence 

fit) = ll^(^)L.<.,rf^ = I' \\Jis,t)t^^^^^ds. 

This equation states that f{t) can be written as the limit of a convex combination of convex 
functions Ui(t) = \\J{si,t)\\i^^^ so / must be convex itself. □ 

In a recent paper (Corollary 8.7 of [21j). the authors prove this property of convexity of the 
geodesic distance in a general setting concerning nonpositively curved symmetric spaces given by a 
quotient of Banach-Lie groups. 

Lemma 3.11. For any a;, y G TCr we have 

dist(e^e'^) - |lln(e-/V^e-/2)||, > \\x - y\\, (9) 

Proof. Take j{t) = e*^, 6{t) — e*** and / as in the previous corollary; we may assume that x,y E HS'\ 
Note that /(O) = 0, hence f(t)/t < /(I) for any < i < 1; hence lim f(t)/t < /(I). Now 

fit)/t= i|lln(e*^/2e-*^e*^/2)||^ = tr([i ln(e*^/2e-*''e*^/2)]2)i^ 

and 

lim - ln(e*^/2g-tygt:./2N ^ d ^.^tx/2^-ty^tx/2^ ^ _ y) = ^ „ □ 

t^o+ t dt 
Corollary 3.12. The inner angles of any geodesic triangle in E add up to at most tt. 
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Proof. Using the invariance of the metric for the action of the group of invertible operators, and 
squaring both sides of inequahty ^ in Lemma 13.111 we obtain the Hyperbohc Cosine Law: 

> + ll-i - 2h+ih-i cos(a,). (10) 

Here li (i=l,2,3) are the sides of any geodesic triangle and is the angle opposite to li. These 
inequalities put together show that one can construct a comparison Euclidean triangle in the affine 
plane with sides k. For this triangle with angles Pi (opposite to the side k) we have If = If^^ + 
— 2/i+i/i„i cos(/3i). This equation together with inequality ([TO)) imply that the angle Pi is bigger 
than ai for i ~ 1,2, 3. Adding the three angles we have ai + 02 + < f3i + f32 + = ■ □ 

Proposition 3.13. The metric space (S, d) is complete with the distance induced by the minimizing 
geodesies. 

Proof. Consider a Cauchy sequence {pn} C S. Again by virtue of inequality ^ of Lemma 13.111 
Xn = ln{pn) is a Cauchy sequence in TYk. Since Hilbert-Schmidt operators are complete with the 
trace norm, there is a vector x G such that x„ ^ a; in the trace norm. Since the inverse map, 
the exponential map, the product and the logarithm arc all analytic maps with respect to the trace 
norm, dist(p„,e^) = || ln(e'^/^e~^''e'^''^)||2 when n — > cxo. □ 

4 Geodesically convex submanifolds 

Definition 4.1. A set M C T, is geodesically convex if for any two given points p,q ^ M, the 
unique geodesic of S joining p to q lies entirely in M . A Riemannian submanifold M C T, is 
complete at p £ M if Exp*^ is defined in the whole tangent space and maps onto Ad. The manifold 
M is complete if it is complete at any point. 

Remark 4.2. The manifold S is complete; moreover, Exp^ is a (analytic) isomorphism of Hr 
with S for each p G S. Other notions of completeness are touchy because, as C. J. Atkin shows 
in [3] and [1], the Hopf-Rinow Theorem does not necessarily hold in infinite dimensional Banach 
manifolds. 

These previous notions are strongly related; it is not hard to see that for any Riemannian 
submanifold M of S, M is geodesically convex if and only if M is complete and totally geodesic. 
On the other hand, it should be clear from the definitions that whenever M is a convex submanifold 
of S, M is nonpositively curved. 

4.1 An intrinsic characterization of convexity 

From now on the term convex stands for the longer geodesically convex. As before [ , ] denotes the 
usual commutator of operators in B(iJ). To deal with convex sets the following definition will be 
useful; assume m C Tim is a real linear space. 

Definition 4.3. We say that m is a Lie triple system if [[a, b],c] £ m for any a,b,c € m. Equiva- 
lently, [x, [x, y]] £ tn whenever x,y € m. 

Note that whenever a, b, c are self-adjoint operators, d = [a, [b, c]] is also a self-adjoint operator. 
So, for any involutive Lie subalgebra of operators a C Tic (in particular: for any associative Banach 
subalgebra), m = 9^e(a) is a Lie triple system in Hw- 
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Assume Af C E is a submanifold such that 1 £ M, and M is geodesic at p = 1. Then TiM is a 
Lie triple system, because the curvature tensor at p = 1 is the restriction to TiM of the curvature 
tensor of E, and 9ii{x, y)z = —\\\x, y], 2;]. In particular, if M is geodesically convex, TiM must be 
a Lie triple system. This weak condition on the tangent space turns out to be strong enough to 
obtain convexity: 

Theorem 4.4. (Mostow-de la Harpe JME^\TTj) Assume m C TYr is a closed subspace, put AI — 
exp(m) C E with the induced topology and Riemannian metric. Assume further that m is a Lie 
triple system. Then for any p,q £ M it holds true that qpq € M . 

Proof. As P. de la Harpe pointed out, the proof of G. D. Mostow for matrices in [53] can be 
translated to Hilbert-Schmidt operators without any modification: we give a sketch of the proof 
here. Assume p = e^, g = e^, and consider the curve e"*-'' = e'^'e^e*^. Then it can be proved that 
a{t) = G{a{t)) with G a Lipschitz map that sends tn into m (this is nontrivial). Since a{Q) = x E m 
and G is a Lipschitz map by the uniqueness of the solutions of ordinary differential equations we 
have a C m. Hence e"'"'^' = qpq G M and the claim follows. □ 

Corollary 4.5. Assume M ~ cxp(m) C E, and m is as in the above theorem. Then M is a closed 
convex submanifold. 

Proof. Take p,q £ M. Then p = e^, g = e^ with x,y £m. If we put r = e~^''^e^e~'^/^, then r G M 
because e~^/^ and e^ are in M. Moreover, z — ln(r) G m. But the unique geodesic of E joining p 
to q is 7(t) = e'^/^e*^e'^/^, hence 7 C M. □ 

Corollary 4.6. Assume m C Tim is a closed, commutative associative Banach subalgebra of Tic. 
Then the manifold M = exp(m) C E is a closed, convex and flat Riemannian submanifold. More- 
over, M is an open subset of m and an abelian Banach-Lie group. 

Proof. The first assertion follows from the fact that m is a Lie triple system. Curvature is given by 
commutators, hence M is flat. Since m is a closed subalgebra, = ^ ^ G m for any a; G m, so 
M C m. That M is open follows from the fact that exp is a isomorphism fCorollarv l3.6p . □ 

If M is flat and geodesic at p = 1, TiAf = m is abelian (by Proposition 13. 4|) . therefore 

Corollary 4.7. Assume M = exp(m) is closed and flat. If M is geodesic at p = 1, then M is a 
convex submanifold. Moreover, M is an abelian Banach-Lie group and an open subset ofm. 

We adopt the usual definition of a symmetric space [16] : 

Definition 4.8. A Riemann-Hilbert manifold M is called a globally symmetric space if each point 
p £ M is an isolated fixed point of an involutive isometry Sp : M —>■ M . The map Sp is called the 
geodesic symmetry. 

Theorem 4.9. Assume M = exp(m) is closed and convex. Then M is a symmetric space; the 
geodesic symmetry at p G M is given by Sp{q) = pq^^p for any q G M . In particular, is a 
symmetric space. 

Proof. Observe that, for p — e^, q — e^, Sp{q) = e^e~^e^; this shows that Sp maps M into M. To 
prove that Sp is an isometry, for any vector w G m consider the geodesic ay of M such that a(0) = q 
and d(0) = v. Then a{t) = ge*?"'" and 

d{sp)q{v) = 4-|t=o(sp o av) = -pq^^vq'^p. 
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Since M has the induced metric, \\pq ^vq ^pW^ _i = by Lemma 12.51 (with g = pq ^). In 

particular, dpSp = —id, so p is an isolated fixed point of Sp for any p G M . □ 

Theorem l4.4l and its corollaries imply that E (as any symmetric space) contains plenty of convex 
sets; in particular 

Remark 4.10. We can embed isometrically any fc-dimensional plane in S as a convex closed 
submanifold: take an orthonormal set of k commuting operators (for instance, fix an orthonormal 
basis {ei}i£M of H and take pi = Ci (E) Ci, i = 1, ■ ■ ■ , k), and consider the exponential of the linear 
span of this set. In the language of symmetric spaces, we are saying that rank (S) = +oo. 

Let I{M) be the group of isometrics of a submanifold M. 

Theorem 4.11. If the submanifold M = exp(m) is closed and convex, then I{M) acts transitively 
on M. 

Proof. Take p = , q = two points in M, v ~ p\n{p^^q) and 7(t) — pe*^ the geodesic joining 
p to q. Note that q = 7(1) = pe^ — e^P p. Consider the curve of isometrics ipt — ^7(4/2) ° 
Then 

ip^[p) = e5''""'e^e-^e^'"^"'e^ = e'"=""'e^ = g. □ 

Remark 4.12. Assume M C S is closed and convex, and 1 G M. Let I{M) be the group of 
isometrics of M . Then, since any isometry ip is uniquely determined by its value at 1 S M and its 
differential dtfi, the set I{M) can be naturally embedded in a Banach space: take ip G I{M) and 
consider 

Note that dJpi is a unitary operator of TiM ~ m (with the natural Hilbert-space structure), so 
there is an inclusion J : I{M) ^ M x U{B{m)) given by the map (p 1-^ {(p(l),dlpi). On the other 
hand, for a given pair {p,u) e Af x U{B{m)), put p>{e^) = p^ exp(u(a;))p2 , (a; e m). It is not hard 
to see that (p is an isometry of M which maps 1 to p, such that dTp-^^ = u . Hence we may identify 
I(M) ~ M X U{B{xn)). 

Remark 4.13. If M — exp(m) is closed and convex, it is geodesic at any p = e^ £ M, so 

TpM = Expp ^(M) = {p^- \n{p-^ qp^^)p^ : q e M} 

(see Remark 12.41) . Since p^ — e^^^ £ m, using Theorem 14.41 we obtain the identification TpM = 
p2 (TiAI)p'^ — p2 tnp2 . It also follows easily that an operator v S Hr is orthogonal to Af at p 
(that is, V e TpM^) if and only if 

<p ^ zp ^ , V >2=< P ^ vp ^ ,z = for any z e m. 

In particular, TiM^ = = {w G TLr : < v,z >2— for any z e m}. Note that, when m is a 
closed commutative associative subalgebra of operators, y p^yp^ is a linear isomorphism of m; 
in this case TpM = m = TiM for any p G A/. This last assertion also follows easily from Corollary 
WM. and clearly TpM^ = TiM^ = in this case. 

Remark 4.14. Assume Af C S is a convex submanifold. If the curve 7 is the geodesic joining p 
to q, then the isometry (pt = 57(4/2) ° Sp translates along 7, namely 

pthis)) = p e^P''" p-'^ p e'P''" p-^ p eiP^'" = 
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In particular, fiip) = q. Now take any tangent vector w € T^(^sjM, and let 

wit) = {dipt)^(^,-){w) = e^^P^'we^P^'". 

It follows from a straightforward computation using equation ^ of Section [5] that w(t) is the 
parallel translation of w from 7(3) to 7(5 + t); namely w = 0. We conclude that the linear 

map {dipt)^(^s) ■ T^{s)M T^i^s+t)M gives parallel translation along 7, i.e {dift)-y{s) = Pl'^'^in)- 

1 - ^ ^1 

particular, since q — ^{1) — p"^ eP ^^p^ p^ , the map 

Pj^ : W 1—^ P^ (p 2qp 2j2p 2 p 2 2gp 2)2p2 

gives parallel translation from TpM to TqM. See also Theorem 14. 181 

4.1.1 Examples of convex sets 

1. For any subspace s C Tin, tUs = {x e : [x, y] = V y G s} is a Lie triple system. 

2. In particular, for any y e Hr, rrij, = {x e TYr : [x, y] = 0} is a Lie triple system. 

3. The family of operators in TYr which act as endomorphisms of a closed subspace S d H form 
a Lie triple system in Hr. 

4. Any norm closed commutative associative subalgebra of TYr, closed under the usual involution 
of operators, is a Lie triple system. In particular 

(a) The diagonal operators (see Section [H]). This is a maximal abehan closed subspace of Wr, 
hence the manifold A (which is the exponential of this set) is a maximal flat submanifold 
of E. 

(b) The scalar manifold A = {Al : A e M>o} C E is the exponential of the Lie triple system 
M 1 C Wr . 

(c) For fixed a € HS'\ the real part of the closed algebra generated by a, which is the closure 
in the 2-norm of the set of polynomials in a, is a Lie triple system. 

5. The real part of any Banach-Lie subalgebra of Tic is a Lie triple system (in particular: the 
real part of any associative Banach subalgebra). 

4.2 Convex manifolds as homogeneous manifolds 

The results of this section are related to those of Sections 3 and 7 of Chapter IV in [TB] . See also 
Theorem 5.5 in [21] for a proof of the existence of smooth polar decompositons in the (broader) 
Banach-Finsler context. 

Definition 4.15. Let 7Y* be the group of invertible elements in Tic- This group has a natural 
structure of manifold as an open set of the associative Banach algebra Tic; it is a Banach-Lie group 
with Banach-Lie algebra Tic- 

LetUiTic) stand for the unitary elements of the involutive Banach algebra TLc, namely the set 
of u E ?i* such that u* — u^^. It is a real Banach-Lie subgroup o/Ti* with Lie algebra iHr. 

Let G be a connected abstract subgroup o/7i* . We say that G is a self-adjoint subgroup of Tic' 
if g* £ G whenever g G G (for short, G* = G). Note that a connected Banach-Lie group G is 
self-adjoint if and only if g* — g, where g denotes the Banach-Lie algebra of G. 
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If a d TLc is a linear space over M, let [a, a] = span^{[a^ 6] : a, 6 S a}, where the bar denotes 
closure in the norm of the Banach algebra Tic- 

If A G 7i* is a set, {A) will denote the abstract subgroup generated by A (the group whose 
elements are the inverses and the finite products of elements in A). 

Let I X 1= {xx*)^ = exp(i ln(a;x*)) for x £ Tic- Since Tic is an involutive Banach algebra, 
\x\eJ:cHlifxe He- 

Remark 4.16. The group He', having the homotopy type of the inductive hmit of the groups 
GL{n,€.) (see jjj, Section II. 6) is connected; moreover, there is a homotopy class equivalence 

He' -S^xS^x SU{oo). 

Here SU{oo) stands for the inductive limit of the groups SU{n,C). 

Proposition 4.17. Let g C Tic be a closed real Banach-Lie subalgebra. Then G = (cxp(g)) admits 
a topology and a smooth structure such that G is a connected real Banach-Lie group and g = TiG 
is the Banach-Lie algebra of G. The inclusion G ^ TiJ is a smooth inmersion and the exponential 
map of G is given by the usual exponential of Tic ■ The topology on G might be strictly finer than 
the topology of Tic- 
Proof. Since Tic is a Hilbert space, the Banach-Lie subalgebra admits a suplement. By Theorem 

j 

5.4 of Chapter VI in [^U], there exists an integral manifold H ^ Ti^ for the subbundle {gBjgeH^- 

The manifold H is connected, and a Banach-Lie group with dji{TiH) = g. Since j is a smooth 
homomorphism of Banach-Lie groups, we have j o Exp^ = cxpodji. The other assertions follow 
from this identity because G = (exp(g)) = o Exp^ (TiH)) — j{II). □ 

Theorem 4.18. Let G — (exp(g)) C Tic* be a connected self-adjoint Banach-Lie group with 
Banach-Lie algebra g C Tic- Let P be the analytic map g i— s- gg*, P . G ^ G. Let { = ker(dPi), 
m = Ran(dFi). Let Mq = exp(m), K = GC\U{Tic) = ^'"^(1). Then 

L The set m is a closed Lie triple system in TYr. We have [tn, m] C i, [m,6] C m, C t and 
g = m©t. In particular, t is a Banach-Lie subalgebra of g (and of iTis. also). 

2. P{G) — Mq, and Mq is a geodesically convex submanifold ofY,. 

3. For any g =\ g \ Ug & G (polar decomposition) , we have \ g |e Mq and Ug G K . 

4. Let g e G, p e Mq, Ig{p) = gpg* - Then Ig e I{Ma)- If g = p^ {p^^qp^^)^p^^ £ G, then 
Ig{p) — q, namely G acts isometrically and transitively on Mq- 

5. Let u G K and x G m (resp. ). Then Iu{x) = uxu* e m (resp. m^). If p,q G Mq then Ip 
maps TqMc (resp. TqMc^ ) isometrically onto Tj^(^g-^MQ (resp. Tj^(^gjMQ'^ ). 

6. The group K is a Banach-Lie subgroup of G with Lie algebra t. 

7. G ~ Mq X K as Hilbert manifolds. In particular K is connected and G/K Mq. 

Proof. 1. Note that dPi{x) — x + x* , hence t — {x G q : x* = —x} which is certainly a closed Lie 
algebra. Note also that m = {a; G g : a;* = a;} is a Lie triple system; it is closed because a; 1— > x* is 
an isometric automorphism of Tic- Since [x,y\ = xy — yx is self-adjoint whenever x is self-adjoint 
and y is skew-adjoint, the other assertions are clear. 

2. Cleary P{G) 13 exp(m) because e^ — P(e^/^). On the other hand, since g splits, there exist 
neighbourhoods of zero Um C m and J/j C 6 such that the map Xm+Vi ^ e^^e^' is an isomorphism 
from Um®Ui onto an open neighbourhood V oil € G. Then {V) is open (and closed) in G and so is 
all of G. Hence, for any g G G, g = (e^^e^^)"^ ■ • ■ (e^"e^")"" for self-adjoint Xi G Um, skew-adjoint 
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Ui € C/(, and = +1. Now e^e^e^ G exp(Tn) whenever a;,?/ G m f Theorem I4.4p . and inspection of 
the expression for P{g) = gg* shows that P{g) lies in exp(m) if e^e^e"^ G exp(m) whenever x G m 
and y G Equivalently, we have to show that ^c?(e^) maps m into m; since Ad{ey) = e°''^^v\ it 
suffices to show that ad{y) : x ^ [y, x] maps m into m, and this follows from the previous assertion. 
The set Mq = exp(m) is a convex submanifold because m is a closed Lie triple system (Corollary 

3. If g G G, then gg* = e"^" for some G m. This implies that | g |= e^°/'^ e Mq C G. Now we 
have Ug =\ g g ^ G, and clearly Ug G K. 

4. If p G Mg, then p = P{go) = gogl for some go G G. Then, if g G G, Ig{p) = ggoglg* = 
P{99o) G Mg- Note that Ig is an isometry of Mg, because Mg has the induced metric, so Lemma 
12.51 apphes. 

5. If a; G m and u ^ K, then G hence ue^u* — exp(J„(e^)) G A/q. Hence uxu* = 
\n{ue^u*) G m. Since < /„(?/), a; >2=< y,Iii'{x) >2 (see Remark [^TT]! . we obtain the proof of the 
assertion concerning rri"^ . 

Clearly Ip maps TqMG isometrically onto Tj^^^q-^MG- Assume now w G TqMG^ — q^m^q^ (see 
Remark I4.13p . If u = {pqp)^p~^^q^^ , then u E G and uu* — 1, hence wq — u{q^^wq^i)u* G 
by the previous assertion. Then /p(w) = pwp — (j)qp)^wo{pqp)2 g Tj^(^q-jMG^ ■ 

6. The previous items show that P : G ^ Mg is surjective. Now dPg : g.Q — > Tgg-MG is given by 
g.x I— > Ig{x + X*). Clearly this map has split kernel gt. Let g —\ g \ Ug a,s above. For z G Tgg*MG 
we have, by Remark I4.13[ z = {gg*)'^w{gg*)^ — I^g^(w) for some it; G m. Let x — (ti;/2), 
then X G m C g and dPg{gx) = z. Hence the group K — P^^(l) is a submanifold of G because 
P : G ^ Mg is a submersion (Proposition 2.3 of Chapter 11 in ■2DJ). 

7. The map T : Afg x K ^ G given by r(p, m) = pu is clearly smooth and it is a bijection by 
the statements above. The inverse is given hy g i-^ {\ g \,\ g g); since | g \— exp(i \n{gg*)), the 
map T is a diffeomorphism. □ 

Remark 4.19. For M — exp(tn) a convex closed manifold in E, consider qm = m © [m, m]. Then 
is a Banach-Lie subalgebra of Tic due to the formal identity 

[[x, y], [z, w]] + [z, [w, [x, y]]] + [w, [[x, y],z]] = 

and the fact that m is a Lie triple system. Let Gm — (exp(gM)}- Then Gj\/ is a connected Banach- 
Lie group with Banach-Lie algebra qm (Proposition 14.17]) . Since (a + [6, c])* = a + [c,b] for any 
a, 6, c G m, then M C Gm and G%f = Gm- It is also clear that 6 = [m, m] {t as in Theorem 14.18^ . 
The elements of M are indeed the positive elements of Gm , and the elements of the stabilizer of 
1 are the unitary operators of Gm- Note that Gm is a submanifold of 7i* if and only if K is a 
submanifold oih({Hc)- 

When m is a commutative associative subalgebra, we have qm = in and also Gm = C m is 
an open set (in particular Gm is a submanifold of HJ). In any case m — Z{m) Z{m)^ = mo ® mi 
(here Z{m) denotes the set {x & m : [x,y] = Vy G m}), and M = exp(mo) exp(mi) ~ AIq x Mi 
where Mi are convex and closed, hence Gm — Mq x Gm^ . Since < x, [a, [b, c]] >2—< c, [b, [a, x]] >2 
for any a, b,c,x E m, it is easy to see that Z{m) = [m, [m, m]]^. 

The results above assert that, for a given convex submanifold M = exp(m), we have Mgm = M. 
On the other hand, for a given connected involutive Banach-Lie subgroup G, we have Gmg C G, 
though in general [m, m] can be strictly smaller than t, so the other inclusion does not necessarily 
hold. The equality holds iff t is semi-simple, i.e. [6, t] — t (equivalently, if Z(£) = 0). 
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It is a well known result (see [H], p.42) that [HS, HS] = HS and [HS'', HS''] = MS''. Therefore 
taking m = HS'', we get 6 = iHS'', and then 0m = HS. This implies Ge/r = Hc'/C^l- Clearly 
-P(Ge) = P{Hc') = S, because any positive invertible operator has an invertible square root. On 
the other hand it is clear that the isotropy group K equals U {He) (the unitary group of He)- So 
there is an analytic isomorphism given by polar decomposition: S 2± Tic* /U (Tic)- The manifold 
of positive invertible operators S is an homogeneous manifold for the group of invertible operators 
Tic* I which acts isometrically and transitively on E. This last statement is well known, and Theorem 
14. 181 can be read as a natural generalization. 



5 Projecting to closed convex submanifolds 

We refer the reader to [2(F for the first and second variation formulas. 

Proposition 5.1. Let M he a convex subset of E, and let p G E. Then there is at most one 
normal geodesic 7 o/ E joining p and M such that £(7) — dist (p, M). In other words, there is at 
most one point q £ M such that dist(p, q) — dist (p, M). 

Proof. Suppose there are two such points, q and r G M, joined by a geodesic 73 G M, such that 
L(7i) — dist(ji,q) — £(72) — dist(p, r) — d{p,M). We construct a proper variation of 7 = 71, 
which we call Tg. The construction follows the figure below, where cr(s,i) := (Js{t) is the geodesic 
joining p with 73(5). 




r(s,t) 



a{s, t) 0<t <1 

73(s(2-<)) l<t<2 



T-.n .^ / 7i(i) if < t < 1 , • / 7i(<) if 0<t<l 

^(i)-mt) = ^^ ^^^^2 . «o that 7(0 = I if ,<,<2 

Also note that the variation vector field (which is a Jacobi field for 7) is given by 



< t < 1 



(2-t)73(0) l<t<2 



If Ai7 denotes the jump of the tangent vector field to 7 at ti, namely jitf) — j{t^ ), and F is a 
proper variation of 7, then the first variation formula for the curve 7 : [0, 2] — > E reads 

II7II. ^|_„+£(r.) = -^' <y(t),A7(i) dt -J2<vit^)^^a>-,iu) ■ 

In this case, Dtj is zero in the whole interval [0, 2], because 7 consists (piecewise) of geodesies. The 
jump points are to = 0, ti = 1 and t2 — 2, so the formula reduces to 



<73(o),7i(i) >,= ^L.o+^(rs)ll7ll,. 
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Recall that 73 C M, and that 71 is minimizing. Then the right hand term is nonnegative, which 
proves that the angle between 71 and 73 at q is bigger that 7r/2. With a similar argument, we 
deduce that the same holds for the angle between 72 and 73 at r. Hence, the sum of the three inner 
angles of this geodesic triangle is at least tt. Since the sum cannot exceed tt (see Corollary I3.12p . 
it follows that the angle subtended at p must be zero, which proves that 71 and 72 are the same 
geodesic, and uniqueness follows. □ 

Next we consider the problem of existence of the minimizing geodesic. 

Proposition 5.2. Let M he a convex submanifold of T,, and p a point of E not in M . Then the 
existence of a geodesic (3 joining p with M such that L{j3) = dist(p, M) is equivalent to the existence 
of a geodesic 7 joining p with M with the property that 7 is orthogonal to M . 

Proof. In fact, the existence of such a geodesic /3 is equivalent to the existence of a point q G M such 
that dist(p, M) — dist(p, q). We will show that if g G M is a point such that jgp is orthogonal to M 
at q, then dist{q,p) — dist(M,p). The other implication follows from the uniqueness theorem above. 
Consider the geodesic triangle generated by p, q and d, where d is any point in M different from q. 
Since jqp is orthogonal to TqM , it is orthogonal to ^qd- Then, by virtue of the Hyperbolic Cosine Law 
(equation ^ in Section^, dist{d,pf = L{jdpf > L{'^qpf + L{-^qdf > L{lqp? = d\st{q,pf. □ 

This last proposition raises the following question: is the normal bundle NM of M diffeomorphic 
to S, via the exponential map? 

Lemma 5.3. Let M be a convex and closed submanifold. Let E : NM T, he the map {q, v) >—>■ 
Exp^{v). For e> 0, put NM^ = {{p, v) £ NM : \\v\\^ < e} and = E{NMe). Then E is injective 
and there exists e > such that E : NM^ — > fi^ is a -diffeomorphism. The set U,^ is an open 
neighbourhood of M in S. 

Proof. Let us prove first that E is injective. Assume there exist p,q G M, v £ TpM-^, w G TqM^ 
with Expp(w) — Exp^(w). Naming r to this point, consider the geodesic triangle in S spanned by 
p,q M and r G S. The geodesic which joins p to r is clearly 71 (t) = Expp(tu), which is orthogonal 
to M at p, and the same is true for 72 (i) = Exp^(tu'), which joins q to r. Hence p — q and v = w 
because of Corollarv l3.12l 

We may assume that 1 e M. Since E{q,v) = qe'> the differential of E at (1,0) £ NM is 
the identity map because TiM © TiM^ = TiE and dexpQ — id. The inverse mapping theorem 
([20], Theorem 5.2 of Chapter I) gives C"-diffeomorphic neighbourhoods Uq — {{q,v) e NM : 
dist(g, 1) < e, ||u||^ < e} C NMe and flo = E{Uo) C S respectively. For given {p,v) £ NM^, 
consider the isometry of M given by /p : x 1— > p^xp'^ , and note that /p(l) — p. If {q,w) G Uo, 
then clearly Lp{q) £ M. Moreover, Lp{w) e Tf^^^^M-^ by Theorem [4.181 hence Up = {Ip x Lp){Uo) 
is an open neighbourhood of {p,v) in NM^ diffeomorphic to Uq. Now E Up — > E{Up) is a 
diffeomorphism, because a straightforward computation shows that E |(7p= LpoEo (Ip x Ip)~^. □ 

Remark 5.4. Clearly E{NM) C S is the set of points p £ S with the following property: there 
is a point q € M such that dist(g,p) — dist {M,p). Note that the map Hm ■ E{NM) M, which 
assigns to p e E{NM) the unique point q G M such that dist(g,p) = dist {M,p), is surjective. This 
map is obtained via a geodesic that joins p and M, and this geodesic is orthogonal to M, therefore 
we call Hjv/ (p) the foot of the perpendicular from M to p. 

Lemma 5.5. Letp, q G E{NM), andIiM{p) 7^ Ha/ (q). If^p is a geodesic that joins Hj\/(p) top and 
7g is a geodesic that joins Hm('?) to q, put f{t) = dist (7p(t), 7g(t)). Then the map f : K>o M>o 
is increasing. 
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Proof. Since / is a convex function (Corollarv l3.10p . it suffices to show that / '(0+) > 0. 

Take a variation a{t,s), where (Tt{s) is 
the geodesic joining jp{t) to jq{t). Then 
a{t,0) = 7p(t), a{t,l) = -fq{t), and 
(t(0, s) — 7(s) is the geodesic joining 
Ha/ {p) to Hm (q) (which is contained in M 
by virtue of the convexity) . Note also that 
(t(1, s) is the geodesic joining p to q. This 
construction is shown in the figure on the 
right. 

Note that f{t) — L{at)- Put V = J^lt^oC- We apply the first variation formula to obtain 

Il7ll. rf^l.^o+^(^t) = -/ <V{s),DM>,^.,ds + 

+ < y(l),7(l) - < ^(0),7(0) >n„(,) . 

The fact that 7 is a geodesic reduces the formula to 

II7II-, /'(0+) = - < ^(1),-7(1) >n,,(.. + < -^(0)'7(0) >n„<„ • 

Note also that V^(0) = 7p(0), V{1) — 7^(0). Recalling that the angles at M are right angles, we 
obtain / '(0+) = 0. □ 

Theorem 5.6. The map Hm is a contraction, namely dist {IlM{p),^M{q)) < dist(p, g). 

Proof. We may assume again that p,q ^ M, and that IIm{p) ^ ^M{q)- In the notation of the 
lemma above, note that /(O) = d {IlM{p),^M{q)) and /(I) = dist(p, q); since / is increasing, the 
assertion is proved. □ 

We want to prove that E{NM) = E. We will do this by proving that it is both open and closed 
in S. The following argument is similar to the one used by H. Porta and L. Recht in 26J. 

Lemma 5.7. For X e [l,+oo), put -qx : E{NM) E{NM), ?7A(Expp(w)) = Expp(Au). Let D,^ be 
as in Lemma \5.3l Then E{NM) = U ri\(Q^), and each rjx : fl^ ^ T, is a diffeomorphism onto 

its open image. 

Proof. Clearly Ux>ir]x{fli) C E{NM). Let us prove the other inclusion. First, if r = Expp(?;) with 
\\v\\^ < e then r £ V,^ — i]i{n^). Let us consider the case where \\v\\^ > e; then r — Expp(w) = 

Expp(Au') with A = ^I^Ip. and w = ^n^, v, so r G r]x{fle) because \\w\\^ — e/2 < e and A > 1. 



Assume that there exist ri,r2 G fie and A > 1 such that rix{ri) — rix{r2)- That is, assume there 
exist p, g e M, V G TpM^, w e TqM^ with \\v\\p < e, ||w||^ < e and Expp(Aw) — Expq{Xw), namely 
E{p, Xv) = E(q, Xw). Since E is injective by Lemma 15.31 we have p — q and v — w. This argument 
proves that the maps r]x are injective. 

Next we show that, for any A > 1 and r G Q^, d (77^)^ : TVS T^^j^.^jS is a linear isomorphism, 
and this will prove the final assertion. Take a C fi^ a geodesic such that q;(0) = r and a{Q) = x. 
Since a is a geodesic, we have that dist(Q;(i), r) = i||d(0)||^ for t > (see Section lXH) . Put /3 — rjxoa. 



Then /3(0) = r^xir) and /3(0) = di7jx)r (x). Clearly dist(/3(t), r?A(r)) < L*(/3) = /J ||/3(s)||^,(^,ds. On 
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the other hand, dist{rix{a{t)),r]\{r)) > dist{a{t),r) — t\\x\\^ where the inequaHty is due to Lemma 
15.51 because A > 1. If we put together these two inequahties and divide by t, we get 

\J^\\mKu,ds>\Mi- 

Taking hmit for t 0+ gives \\d{T]x)r. i^)\\,,^(r) ^ II^L- Now put Ax = /^^^^^^ °d{i]x)^ o Ir, where 
Ip : V t-^ p'^vp'^ are hnear isomorphisms (see Lemma l2.5p . If we consider Ax '■ TiE TiT, = 
the last inequahty says that ||AA(a;)||2 > \\x\\^ for any x G Hr. 

Clearly 771 = id^^ and rf(?7i)^ = idx^^- Since the map (A, r) fjxir) is analytic from R>o x fl^ 
to S, there is an open neighbourhood of 1 G M such that Ax is an isomorphism. Assume Ax is 
invertible for A G [l,m): then 11^4^^ ||b(-h,) < 1 for any A G [l,m). Since Am = lim Ax (in the 

A — >7n^ 

operator norm of B{Hr)) and — 1|| < \\Am — < 1 if A is close enough to m, it follows 

that AmA^^ is invertible, thus Am is invertible. Since the maps Ip are isomorphisms, d{rjx)^ is an 
isomorphism for any A > 1, and any r G ri^. □ 

Corollary 5.8. The set E{NM) is open in S. 

Theorem 5.9. Let M be a convex closed suhmanifold of S. Then for every point p G S, there 
is a unique normal geodesic jp joining p to M such that L(^p) = dist(p, M). This geodesic is 
orthogonal to M, and ifllM : S ^ M is the map that assigns to p the end-point ofjp, then Hm is 
a contraction for the geodesic distance. 

Proof. The theorem will follow if we prove that E{NM) — E. Since E is connected and E(NM) 
is open, it suffices to prove that E{NM) is also closed. Let p G E{NM). There exist points 
Qn G M, Vn G Tq^M^ such that p — limpn — lim Exp (w„). Observe that g„ — IiM{Pn), 
SO dist(g„,g„i) < dist(p„,Pm)- Since {pn} converges to p, it is a Cauchy sequence. It follows 
that {qn} is also a Cauchy sequence. Since M is closed (and therefore complete), there exists 
q E M such that q = limg„. We assert that dist(p, g) ~ dist(p, M). First note that dist(p, <7„) < 

n 

dist(p,p„) +dist(p„,g„) and dist(p„,g„) = dist(p„,M), so dist(p, g„) < dist(p,p„) +dist(p„,M). 
Taking limits gives dist(p, q) < dist(p, M). □ 

Note that E decomposes as a direct product: with the contraction IIm, we can decompose E by 
picking, for fixed p, 

1. the unique point q — IlAf (p) such that dist(p, q) — dist(p, M) 

2. a vector Vp normal to TqM such that the geodesic in E with initial velocity Vp starting at q 
passes through p; note that Vp ~ Expj^^^(.p^ (p), and also \\vp\\^ — dist(p, M). 

Since the exponential map is analytic on both of its variables, we get 

Theorem 5.10. The map p ^ (IlMip), Vp) is the inverse of the map (g, Vq) Exp^(ug), and gives 
a real-analytic isomorphism between the manifolds E and NM . 

Theorem 5.11. Fix a closed convex submanifold M ofT,. Let a G E. Then there exist unique 
operators c G E, w G such that c G M, v G TcM-^ , and a — ce'^ 

Using the tools of Section |4l we can write the factorization theorem in terms of intrinsic operator 
equations (see p3^ for the finite dimensional analogue): 
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Theorem 5.12. Assume m C TYr is a Lie triple system. Then for any operator a G TYr, there exist 
unique operators x € m and v G xn'^ such that the following decomposition holds: e" = e" e^. The 
map da '■ TYr — > R, da{y) = \\ ln(e''/^e~^e''/^) has the operator 2x as its unique minimizer in m. 

As a corollary, we obtain a polar decomposition relative to a convex submanifold. 

Theorem 5.13. Assume AI — exp(m) C S is a closed convex submanifold. Then for any g £ Tic* 

there is a unique factorization of the form g — e^e^u where x G m, v € and u S U{'Hc) is a 
unitary operator. The map g i— » (6^,6", m) is an analytic bijection which gives an isomorphism 

He' ~ X exp(m^) x U{nc). 

Proof. Since gg* G E, we can write gg* = e^e^^e^ with x £ m and v e m^. If u = {e^e")~^g = 
Q-fQ-Xg have uu* — e""" gg* e""" — 1 and also u*u — g*e~^e~^e^^e^^g ~ 1. Hence u is a 
unitary operator and g = e^e^u. This factorization is unique because if g = e^^e^^wi = e^^e"^U2, 
then gg* = e'^^e^^^e^^ = e^^e^^^e^^, so xi = X2, vi = V2 and then m — U2. □ 

6 Projecting to the manifold of diagonal operators 

Lemma 6.1. Let a, /3 G M and a,b £ HS^ . Then 

Exp„+,(/3 + 6) = ae/'/" + fc 

where k is a self-adjoint Hilbert- Schmidt operator. 
Proof. It is a straightforward computation: 

(a + a)e("+°)"'(''+^) = (a + a)[l + (a + a)^i(/3 + b) + ■ ■ ■] 

= {a + a)[l + /3/a+^{P/af + --- + k]. □ 

We need some remarks before we proceed. Fix an orthonormal basis {cijigN of H. 

1. Consider the diagonal manifold A C S: 

A = {d + a > : a £ S., rfisa diagonal Hilbert-Schmidt operator}. 

It is closed and geodesically convex. This is due to the fact that the diagonal operators form 
a closed commutative associative subalgebra. 

2. If do G A, then Td„A = {a + d; a E M., d G HS is diagonal and real} = TiA (see Remark 
EM- 

3. Consider the map A i-^ A^ — the diagonal part of A. Then 

(a) For Hilbert-Schmidt operators we have A^ — J^iPi^Pi where convergence is in the 2- 
norm (and hence in the operator norm); here pi = =< e^, • > is the orthogonal 
projection onto the real line generated by e^. 

(b) (A^)^ = A^ and tr {A'^A) = tr{{A^)^). 

(c) tr{A^B) = tr{AB) if B is diagonal. 

4. The scalar manifold A = {Al : A G K>o} is convex and closed in S, with tangent space at any 
A G A given by Ml C Hr. 
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5. A vector u = /i + u is contained in T(JqA^ if and only if /i = and = 0. This follows from 
Remark 14.131 the fact that + G Td^A, and Remark (3) of this list. In other words for 
any do G A, 

Trf„A^ = TiA^ = {v e HS'' : v is codiagonal } T. 
Theorem 6.2. Let a £ HS'\ Then there exist A £ M>o, d £ A and a; £ HS'* such that: 

a + A = (d + A)e('*+^)~'^ = (d + A)^e('^+^)"^"('*+^)"^ (d + A)i 

Moreover, for fixed \, d and v are unique and a + X ^ id,v) (which maps E — > NA) is a real 
analytic isomorphism between manifolds. 

Proof. Let A = ||a||oo + e, for any e > 0. Then p = a + A £ E. Let IIa{p) = d + a, where d £ A. 
Now pick the unique v £ Td+aA^ such that Exp^^_^(u) = p, this operator v has the desired form 
because of Remark (5) above. As a consequence of Lemma I5TT] a = A, for in this case (3 = 0. □ 

This theorem can be rephrased saying that, given a self-adjoint Hilbert-Schmidt operator a, 
for any A £ ]R>o such that a + A > 0, one has a unique factorization a + A = D e^D where 
D = {X + d)i > is a diagonal operator and w = D~^vD~^ € F is a self-adjoint operator with null 
diagonal. The normal bundle clearly splits in this case, so 

Proposition 6.3. Consider the submanifolds A, exp(r) C S. Then the projection map Ha induces 
a difjeomorphism S ~ A x exp(r). 

Corollary 6.4. For any g £ Tic*; there is a unique factorization g = de^u, where d is a positive 
invertible diagonal operator of Tic, w is a self-adjoint operator with null diagonal in Tic and u is a 
unitary operator of Tic ■ 

Proof. The previous results together with Theorem 15. 131 □ 

7 A foliation of codimension one 

In this section we describe a foliation of the total manifold, and show how to translate the results 
from previous sections to a particular leaf (the submanifold Ei) in order to show an aplication 
concerning (finite dimensional) matrix algebras. Recall that we write HS'' for the self-adjoint 
Hilbert-Schmidt operators. Fix A £ R>o. Let 

Ea = {a + A £ E, a £ HS''}. 

Observe that E^ H = when A 7^ /3, since a + X = b + (3 implies a — b ~ (3 ~ X. In this way, we 
can decompose the total space by means of these leaves, E = U Ea- 

Proposition 7.1. The leaves Ea are geodesically convex closed submanifolds. 

Proof. We consider the projection to the convex scalar manifold A (see Remark (4) above). The 
fact that the projection IIa is a contraction (therefore a continuous map) implies that Ea is closed; 
one must only observe that Ea = 11^^ (A). To show that Ea is geodesically convex we recall that, 
by virtue of Lemma l6.ll for any real A > and any p £ Ea, there is an identification via the inverse 
exponential map at p, TpY}\ = HS'*. □ 
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Remark 7.2. Take (5 + c G Ta+x'Sx^. Since Ta+\'E\ can be identified with HS'', the equahty 
< 5 + c,d V d e HS'' is equivalent to 

tr [{a + X)-^ [{6 + c){a + Xy^ - 5/X] d] ^ V d e HS''. 

Equivalently, T^+a^a ~ span{a + A); shortly TpY.^^ = span{p) for any p G Sa- 

Proposition 7.3. Fix real a, A > 0. Let Ila,\ = He^ |s„ : Sa ^ Sa- T/ien 

1. nQ,A(p) = ^P, so Ha.xip) commutes with p. 

2. Ila,x is an isometric bijection between and Ea, with inverse IlA^a. 

3- He A gives parallel translation along vertical geodesies joining both leaves (that is, geodesies 
orthogonal to both leaves). 

Proof. Notice that for a point b+a G Sq, to be the endpoint of the geodesic 7, starting at a+X G Ea) 
such that £(7) = dist (b + a, Sa), we must have 

b + a = Exp„+^(a; + c) = Exp^^^{k.{a + A)) = e'=(a + A) 

where G R comes from Remark 17.21 above, since x + c G Ta+A^A^- From Lemma [6. 1[ we deduce 
that fc = In ( j), and a = ^b. So, b + a = j(a + A) and also j{t) — (a + X) (j) . Now it is clear 
that Ha (6 + a) — ^{b + a) commutes with b + a. To prove that 11 is isometric, observe that 

dist(n„,A(p),n„,A(<?)) = \\Hi-prH-q){-pr'')\U = ||in(p-^gp-^)||, =dist(p,g). 

a a a 

That n gives parallel translation along 7 follows from the formula for 11 given in the first item of 
this proposition and the formula for the parallel translation given in Remark 14. 141 □ 

The normal bundle in the case of M — Ei can be thought of as a direct product: 

Proposition 7.4. The map T : E ^ Ei x A , which assigns a + a >— > (^(a + a), a) is bijective and 
isometric (T,i and A have the induced submanifold metric). In other words, there is a Riemannian 
isomorphism E ~ Ei x A. 

Proposition 7.5. The leaves E^, Ea are also parallel in the following sense: any minimizing 
geodesic joining a point in one of them with its projection in the other is orthogonal to both of them. 
See Figure 1 below. 



5(t)=c'-l/S^ 



A 



/ • 


V 










1 r 


T 



L(y)=L(S)=dist(Z..XJ 



Figure 1: The geodesies 7 and S are minimizing, the geodesic /3 is not. 
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For any 6 + a G we have dist(& + a, S^) = dist(SQ,,EA) =| lii(j) |. In particular, the 
distance between a, A in the scalar manifold A is given by the Haar measure of the open interval 
(a,/3) on the multiplicative group M>o- 

Proof. It is a straightforward computation that follows from the previous results; the last statement 
was observed by E. Vesentini in another context [28]. □ 

Since E is a symmetric space, curvature is preserved when we parallel-translate bidimensional 
planes. Note also that vertical planes (i.e. planes generated by a vector v e HS'' = T\T.x and A) 
are commuting sets of operators. 

Proposition 7.6. Let p G E^. Then the sectional curvature of vertical 2-planes is zero. 

Proof. It follows from the formula for the curvature given in Section 13.11 □ 

7.1 The embedding of M+ in Si 

Let M+ be the set of positive invertiblc n x n matrices (see the introduction of this paper). First 
note that we can embed ^ Ei for any n G N: fix an orthonormal basis {e„}„gN of H, let 
Pij — Ci® Cj , and identify the set AI„ of real n x n matrices with the set 

n 

We identify the manifold Af+ with V — {e^ : T G T} C Ei and the tangent space at each ^ V 
is T. The set V is closed and convex in Ei by Corollary 14.51 Let us call S = span(ei, • • • , e„), 
S-^ — span(e„+i, e„_|_2 • • •)■ The operator P^ is the orthogonal projection to S and = 1 — P^ is 
the orthogonal projection to S-^. Using matrix blocks, for any operator A G 6(5*), we identify 

Remark 7.7. There is a direct sum decomposition of HS'' ~ T (B J' where operators in J E J 
are such that JPg = 0. A straightforward computation using the matrix-block representation 
shows that < a, 6 >2= for any a G T, & G J7, which says T"*- = J (here we consider HS'' as the 
total space). So the manifolds exp(J') and V — exp(T) are orthogonal at 1, the unique intersection 
point. In the notation of Theorem 14. 181 it is also clear that V — Mq ~ G/K, where 

Theorem 7.8. Let V ~ A/+ C Ei with the above identification. Then for any positive invertible 
operator e** G Ei, (b G HS''J there is a unique factorization of the form 

»-(r x)} /-(o s)- 

then = e^P^ eP^ M+ , X* = X is a Hilbert- Schmidt opertor acting on the Hilbert space 

S^, andY e B{S,S^). 
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An equivalent expression for the factorization is 

b ( e^/' \ r / e-^/2y* \ 1 / e^/2 

ij^^PlUe-^/^ X j/U 1 

Proof. From previous theorems and the observations we made, we know that e*" = e^/^Ce"/^, where 

for some A S 6(5*) and some v G HS''. That Vii ^0 follows from the fact ( see Remark 1 7. 7p that 
r-L = J, and V e TeaT'-L if and only if trie-'^Bc-^Vii) = for any B €T. □ 

Remark 7.9. For any b G HS'', the operator 



is the 'first block' n x n matrix which is closest to e*" in E, and with a slight abuse of notation for 
the traces of B(5) and B(5^), we have 



distlP.e") = dist(e^e^) = ^\\Ye-^/X + 

Corollary 7.10. For any g £ Tic' there is a unique factorization g — \re'"u, where A G K>o, 
u G U{TCc) is a unitary operator, 

R \ /or* 
1 y \Y X 

with R G B{S)+ ~ M+, X = X* e B{S^) a Hilhert Schmidt, and Y G B(5, S^). 

Proof. We use the notation of Remark 1 7. 71 Note that, by Theorem 1 5. 131 g = re^u with u G U{7ic), 
r E V = exp(T) and x G . But = J' Ml if we consider TYr as the total space, and 

ga+a ^ gaga if q, g □ 
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